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ABSTRACT 



Pease (1987) related the effects of atmospheric forcing, mainly temperature and wind 
fields, to the size of coastal polynas. Using Pease's formulation and Kilhvorth's (1977) 
plume model as applied by Melling and Lewis (19S2), the effects of atmosoheric forcing 
on brine injection into the Arctic pycnoeline are investigated. This paper will discuss the 
likelihood of coastal polynas as a source for -denser abyssal waters. 

A standard case was developed for the model with initial conditions taken from 
Melling and Lewis (19S2) and Pease (1987) for comparison with individual sensitivity 
experiments. Ten environmental parameters were individually examined for their influ- 
ence on the plume depth after 90 days. The standard case resulted in a 90-day plume 
depth of 436 meters. A submarine canyon case was simulated, resulting in plume pene- 
tration to over 1300 meters in 90 days. Further experiments used actual T-S soundings 
from Aagaard et al. (19S1) and Ostlund et al. (1987). Finally, a 20 kilometer wide plume 
is shown to penetrate to almost 600 meters in 90 days. 
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I. INTRODUCTION 



The Arctic Ocean, a small and remotely located part of the world's ocean, has de- 
manded the attention of those who wish to understand the Earth's climate. The renewed 
emphasis in recent years given to understanding the Earth's climate has underscored the 
role of the Arctic Ocean in the scheme of the global energy balance. The Arctic Ocean 
is also a key player in the world oceans' thermohaline circulation. Because of the im- 
portance of the Arctic Ocean in the global picture, understanding the nuances of the 
Arctic Ocean's physical oceanography and its sensitivities to environmental change is 
essential before climatologists, oceanographers, and meteorologists can develop a more 
realistic and longer range global climatological concept. Two puzzling aspects of the 
Arctic Ocean are its isothermal halocline and the more saline deep waters. How is the 
halocline maintained in a nearly horizontally independent state? What is the mechanism 
behind the formation of the Arctic Ocean's deep waters? These topics will be discussed 
in this thesis. 

Throughout the Arctic Ocean, the characteristics of the upper layers are remarkably 
independent of the horizontal location of the sample. The water column consists of a 
mixed layer, halocline, thermocline, and deeper waters. The mixed layer is a relatively 
fresh water layer approximately 30 meters deep. River run-off and the summer cycle of 
melting ice are considered the main contributors to the fresh water layer. The halocline 
is the layer of water beneath the mixed layer. In this layer, salinities rise from 31 to as 
high as 34.6 ppt at the 200 meter depth level ( Kill worth and Smith, 1984). The halocline 
is. in essence, an isothermal layer with temperatures just above the freezing point. The 
thermocline is a layer below the halocline. in which temperatures rise to approximately 
0° C at 250 meters (Melling and Lewis. 1982) and continues to rise to 0.5 s C at 400 meters 
due to the layer of warm Atlantic water. Below this depth, temperatures decrease with 
increasing depth. The deep waters of the Arctic Ocean are among the most dense waters 
found in the world's oceans and have been suggested as being an additional source for 
ventilating the world's oceans' deep thermohaline circulation (Aagaard,1981). Although 
the physical oceanography of the Arctic Ocean is well described in literature, the exact 
mechanism(s) for the maintenance of the halocline and deep water formation are not 
fully understood. 
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The halocline is generally isothermal except in the Canadian Basin, where warm 
Pacific waters entering via the Bering Strait causes a warm layer at about the 75-100 
meter depths, and presents a paradox when examining its origin. Sandwiched between 
a cold, relatively fresh layer and a warm, saline layer, an obvious solution to the main- 
tenance of the halocline might be simple vertical mixing between the two layers. How- 
ever, vertical mixing of the cold, less saline surface water with the warmer, more saline 
Atlantic water is not considered to be a source since heat is not mixed along with salt. 
Aagaard et al. (1981) showed a discontinuity existed in the T-S correlation slope at 
about 150 meters. Either salinization of surface waters by freezing (i.e. brine rejection) 
or cooling of Atlantic water could be two possible mechanisms to produce halocline 
water. The authors (Aagaard. et al., 1 9 S 1 ) demonstrate mixing between the two layers 
cannot be a source of halocline water by examining the required salinity of the surface 
water. Using a T-S diagram, and assuming linear mixing at a ratio of 3:2 (surface water 
to Atlantic water), the surface water would need to have a salinity of 34.7 ppt. Salinity 
this high has not been observed in the Arctic surface layer. Therefore, the maintenance 
of the Arctic halocline is generally thought to be attributable to lateral advection of cold, 
saline water from the continental shelves. The shelf waters are made more saline due to 
brine rejection associated with ice growth. Aagaard (1981) proposed that continental 
shelf areas where an annual growth of ice less than 2 meters is necessary to modify 
summertime salinity profiles to wintertime profiles may be the source areas of the 
halocline waters. Upwelling of Atlantic water onto the shelves and subsequent cooling 
has also been suggested as a possible source. Upwelling along submarine canyons has 
been reported as early as 1929 (Sverdrup, 1929). As further evidence eluding to the 
halocline's possible shelfwater origin, water on the continental shelves has been shown 
to have temperature and salinity characteristics similar to the waters found offshore at 
a depth of 150 meters (Melling and Lewis, 1982). 

In their paper, Melling and Lewis (19S2) discuss the lateral advection of low tem- 
perature high salinity shelfwaters into the halocline, with particular attention paid to the 
effects of brine rejection due to the growth of ice. Aagaard (1981) used two methods 
of computing the necessary production rate of high salinity shelfwater to maintain the 
halocline. By considering the production rate required to renew the halocline in 10 years 
and necessary mixing ratios to achieve desired T-S characteristics, Aagaard suggests a 
2.5 .y 10 6 to 4.0.vl0-m 3 — sec -1 for a production rate of cold saline shelfwater, assuming a 
halocline renewal period of 10 years. In order to increase the ambient shelfwater salinity 



to values which would form sufficiently dense water to descend the slope into the 
halocline, Melling and Lewis (1982) emphasize the important role of ice dynamics, spe- 
cifically the role of divergence. Areas of large wintertime ice cover divergence allow 
rapid ice growth hence increased salinization of the underlying water column. Leads and 
polynas are locations of rapid ice growth due to the divergence of the insulating ice 
cover. 

Schumacher et al. ( 1983) examined the degree of salinization of the underlying water 
column beneath a polyna. In their study, measurements of salinity and water current 
velocity beneath a polyna located at St. Lawrence Island indicated polynas to be an area 
of high ice growth rates resulting in density driven currents due to increased salinization 
(Schumacher, et al.. 1983). Coupling the salt source, i.c. polyna, with the interleaving 
plumes of dense shelfwater found by Melling and Lewis (1982) may provide a possible 
mechanism for the maintenance of the Arctic halocline. The polyna provides an atmo- 
spheric dependent source of salt into the water column possibly forming plumes which 
descend down into the halocline, as described by Melling and Lewis (19S2). Along with 
the mechanism associated with the halocline's maintenance, the formation of the 
Canadian Basin's deep water has also been attributed to having a source of low tem- 
perature, high salinity shelfwater. 

The dense waters which ventilate the world’s oceans flow out of the Arctic Ocean 
over the Greenland-Scotland ridge system at a depth of about 2600 meters. A single 
source was advocated by Helland-Hansen and Nansen (1909) in which deep water is 
formed in the Greenland Sea. Nansen ( 1906) had written earlier of the possibility of the 
Greenland Sea Deep Water being modified by cooling in the Barents Sea. This theory 
omitted the possibility of deep water formation within the two basins of the Arctic 
Ocean, namely the Eurasian and Canadian Basins. In 1981, Aagaard postulated that 
deep water is also formed within the Arctic Ocean due to the salinity distribution. 
Greenland Sea Deep Water has a salinity of approximately 34.90 ppt. which is the least 
saline of Arctic deep waters. Therefore, at least two sources of deep water formation 
are thought to exist. One source is the traditional Greenland Sea formation area and the 
new source proposed by Aagaard, is within the Arctic Ocean's two basins. 

The characteristics of water mass properties found in the Canadian and Eurasian 
Basins are similar until about the depth of the Lomonosov Ridge. At about 1500 meters 
depth, the Canadian Basin becomes more saline and is slightly warmer than the Eurasian 
Basin Deep Water (Aagaard. 1981). Under the traditional viewpoint of Greenland Sea 
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Deep Water Tilling the two basins through lateral advection, salinities should not exceed 
34.9 ppt. However, during the LOREX (Lomonosov Ridge [Experiment; Weber, 1979) 
salinities in the b . : ,ns vvere determined to be greater than 34.9 ppt. This suggests an 
additional source of salt if not a completely new source of deep water. Current meters, 
placed o' \'r the Lomonosov Ridge, gave evidence to the existence of some tidal and 
other episodic surges carrying salini/cd Greenland Sea Deep Water over the ridge into 
the Canadian Basin. Even with the surges over the Lomonosov Ridge into the Canadian 
Basin, water from the Eurasian Basin must have an additional salt source to obtain the 
characteristics of water found in the Canadian Basin (Aagaard, 1981). Residence times 
can be calculated by comparing the mixing ratios of shelf water to Atlantic Intermediate 
water necessary to achieve the T-S characteristics found in the Canadian Basin Deep 
Water. Considering a small flux of shelfwater (.003 to .06 Sv) necessary to ventilate the 
Canadian Basin over several hundred years may also explain the warmer temperatures 
found in the Canadian Basin (Aagaard ct ah, 1985). By using tracers such as (5 lf, 0 and 
l4 C, Ostlund et al. (1987) agree with the magnitudes of the renewel time scales for both 
basins presented by Aagaard et al.( 1 985). Ostlund et al. (1987) and Aagaard et al. (1985) 
place the residence time of the Eurasian Basin in the order of decades vice several hun- 
dred years (700 years, Ostlund, ct ah, 1987) for the Canadian Basin. The requirement 
remains of an additional salt source at small fluxes to ventilate the Canadian Basin over 
several hundred years point to an episodic mechanism. Killworth and Smith (1984) write 
that shelfwater sinking at the basin edges may occur but a source of waters sufficiently 
dense to penetrate the Arctic Ocean's stable mid-layer remains a problem. Again, using 
a coastal polyna as a salt source, could shelfwater be modified under extreme conditions 
to produce such dense water which would then descend to the deeper depths via gravity 
flow plumes? 

This thesis will combine Pease's numerical formulations relating atmospheric forc- 
ing, mainly wind and air temperature, to the opening and closing of coastal polynas, with 
a streamtube gravity flow model developed by Smith (1975). With simplified 
hydrography and suitable environmental conditions, this model can produce events of 
deep plume penetration along with a more consistent source of high salinity, low 
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temperature water for the maintenance of the Arctic Ocean's halocline. Chapter 2 will 
describe the numerical model used in this study along with an examination into 
drawbacks associated with Smith's streamtube model as described by other studies, i.e. 
Smith (1975). Killworth (1977). and Melling and Lewis (19S2). Chapter 3 will describe 
results obtained by altering environmental inputs into the model and Chapter 4 will 
discuss the results and implications. 



II. NUMERICAL METHODS AND PARAMETERS 



A. THE NUMERICAL MODEL 

With the increased desire to understand the nature of the Arctic halocline and its 
origin, there existed sufficient pieces of the puzzle to develop a coupled model. The 3-D 
turbulent entraining plume model developed by Smith (1975) has been used by several 
studies (Killworth ( 1977). Kilhvorth and Carmack ( 197S). and Melling and Lewis ( 1982)) 
to gain insight into gravity flow currents arising from brine rejection into rhe water col- 
umn due to ice growth. Melling and Lewis (19S2), compute the initial density of the 
plume by estimating the average ice grown annually in a shallow shelf and and the re- 
sultant salinity contribution. In this study. Pease's (19S7) formulation relating a 
polyna's size to offshore wind speed and air temperature, provides a specific source of 
brine which is then coupled with Smith's plume model, as used by Melling and Lewis 
( 19S2), in order to study the possible contribution of the polyna produced denser waters 
to the Arctic halocline and possible formation of Canadian Basin Deep Water. 

1. The coastal polyna 

Pease (1987) describes the polyna's maximum width as an equilibrium between 
frazil ice production and its advection rate. The frazil ice production rate is a function 
of the heat exchange at the air-sea interface. Therefore, with wind speed and air tem- 
perature as independent variables, one can compute the size of the polyna and its frazil 
ice production rate. 

The polyna size can be found by solving: 



(IX 

dt 



= Vi~ 



XpF, 

II; 



( 1 ) 



where is the change in polyna width with time. V, is the advection rate of frazil ice, 
X. is the polyna width, F, is the frazil ice production rate, and H, is the collection depth 
of grease ice or nilas. The solution of equation (1) assumes the meteorological condi- 
tions are steady throughout the polyna's opening time, thus allowing treatment as a 
linear differential equation. Furthermore, if the meteorological event is assumed to be 
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of sufficient duration so that the polyna reaches maximum width, equation (1) can be 
written as: 



, I 11 ; 

X .{maximum ) = — j— . (2) 

* i 

This allows solving equation (2) fox the maximum width after computing the 
frazil ice production rate. F r In this study, meteorological events of interest are 4-10 
days duration, therefore this assumption allowing equation (2) is valid. Ou (19SS) de- 
scribes a more detailed coastal polyna model and concludes the ice edge is less affected 
by higher frequency atmospheric disturbances than longer period, synoptic type vari- 
ations. Although Ou (1988) includes additional physics to model the temporal changes 
of the polyna's ice edge, the steady state solution is still that as given by Pease (1987) 
and shown above in equation (2). 

The frazil ice production rate is given by Pease (1987) as: 



l r7e gTg ~ Qlu + Pa£h£pK^ Tg ~ 

PiL 



(3) 



where the evaporative heat flux has been neglected due to its small contribution relative 
to the uncertainty of the sensible heat flux. Q, u is the upward longwave radiation, aej^ 
is the downward longwave radiation, L is the latent heat of freezing for salt water, p, is 
the density of young sea ice. and p p C h C p V£T a — T„) is the sensible heat flux. 

2. Salinity distribution within the water column 

The salinity contribution to the water column from frazil ice growth is assumed 
to be uniform with depth and is computed using: 



— (f> — S t )F,A. 

S = S' + , — 1 ' 

^new ° < 



{ qdz) 



(4) 
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from Killworth and Smith (19S4). S rew is the new salinity value of the water column after 
brine rejection. 5 is the average salinity of the shelfwater, 5. is the salinity of the ice, A. 
is the effective area of salinization, and f qdz is the areal flux of the underlying water. 
Using F, from the polyna model, one can compute the salinity of shelfwater exiting a 
polyna. 

3. Equation of state 

The density of the current is computed using the International one atmosphere 
equation of state developed by Millero and Poisson (1981). The equation is given below. 

p = p 0 4- AS 4- BS 4- CS (5) 



where, 

A = 8.24493 x 10"’ - 4.0899 x 10 _3 r 4- 7.6438 x 10 _5 r 2 - 8.2467 x 10" V + 5.3875 x 10'V 
B = -5.72466 x 10“ 3 4- 1.0227 x 10~ 4 r - 1.6546 x 10 - 6r 2 
C = 4.83 14 x 10' 4 

and p 0 is the standard density of seawater. This allows the calculation of the current's 
relative buoyancy using: 



A. 



A Po ’ 



( 9 ) 



In this fashion, the results from the polyna model can be used as inputs to the 
streamtube model. 

4. The streamtube model 

Smith's (1975) model of boundary currents includes the effects of friction, 
entrainment, and Coriolis forces. The 3-D plume was modified by Killworth (1977) by 
altering the entrainment rate. Killworth (1977) related the entrainment rate to the 
product of the plume's area and velocity vice only the velocity (Smith, 1975). Melling 
and Lewis ( 19S2) used Killworth's version of the plume model with an entrainment co- 
efficient dependent on the Richardson number. The model used in this paper is the 
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Melling and Lewis (1982) adaptation with the inclusion of several stratified layers to al- 
low a single model initialization. 

Smith's plume model uses two coordinate systems. The Cartesian system is 
aligned with the bottom so the X-axis lies along the shore and the Y-axis lies along the 
slope. The second coordinate system is the curvilinear system with 2 and ;/ axes. The 
position of the plume axis can be defined by the value of c . The angle which the plume 
crosses the isobaths is given by the angle /? . The coordinate system is illustrated in 
Figure 1 on page 10 and includes gravitational, g, and rotational, Q, vectors. 

The equations for the path flow aref 



cent to the bottom. Pressure gradient forces pull the plume down the slope while 
Coriolis forces bend the flow to the right. If frictional forces were not considered, the 
flow would move horizontally along the isobaths (Killworth. 1977) in geostrophic bal- 
ance. Friction is included, using quadratic drag laws to relate frictional resistance to the 
square of the mean velocity (Smith, 1975). The plume is thus allowed to flow down the 
slope, constantly entraining ambient fluid and losing its negative buoyancy. Once the 
plume reaches neutral buoyancy, the plume interleaves with the ambient waters. The 
equations of the plume's flow from Killworth (1977), are: 




( 10 ) 



and 




( 11 ) 



The plume is considered to be a well-mixed, broad, thin layer of high density fluid adja- 



■^r(An = E ] Ai: V, 

“S 




(13) 



-jT- (A VA) = - A VN 1 sin Q sin /?, 



(14) 
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A sin 9 cos P — JV. 



(15) 



These equations are, the along-stream derivatives of the mass flux (12), the momentum 
flux (13), the buoyancy flux (14) and the cross stream momentum balance (15). In the 
above equations, A is the cross-sectional area of the plume, V is the plume's velocity, A 
is the plume s buoyancy, 9 is the slope angle, p is the angle at which the plume crosses 
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the isobath, N is the Brunt- Viiisala frequency, K x is the drag coefficient, /is the Coriolis 
parameter, and £, is the entrainment coefficient. 

The equations were integrated using a modified Adams method. 

Predictor: 



Corrector: 



37 + 1 = >7 + "Tf (^/ - 59/_, + 37 'fi_ 2 - 9/_ 3 ) 



(16) 



37+i = >V + “ 24 (9 / +1 + 19/- o/_, +/_ 2 ), 



(16a) 



with the initial steps provided with the Runge-Kutta method, 

>7+i =>'i + y (*t + 4/S: 2 +^ 3 ) 

*1 = A - X ,■•})) 

^2 =A x i + y • >7 + y ^1) 

^3 =/(*,- + ^->v ~ m») 



(IS) 



This model uses Killworth's (1977) entrainment and frictional 

parameterizations, which has the entrainment rate and frictional force porportional to 
the product of the area and velocity (velocity squared for friction). Thus, 



E 

= ( — )t , where w is the width of the plume and b is the depth of the plume. 
E\ 



K 



w , 1 



Also used in this model is Bo Pedersen's (19S0) entrainment parameterization for small 
slopes. Melling and Lewis (19S2) included this in their version of Smith's model. The 

E = 0.072 sin 6 sin p = 0.072) -§r ) where 

Ki 



entrainment coefficient is given as 
Ri = bA cos d/V 2 is the bulk Richardson number of the How. 

5. The coupled model 

The model takes the form illustrated in Figure 2 on page 13. Underneath the 
polyna, a plume with a depth, b. and a width, w, is salinized due to brine rejection. The 
plume has an initial velocity, V„ . The effects of varying the initial velocity will be ex- 
amined in Section 3, however, the standard case value of V„ = .OAm/s is used. This value 
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is dose to the values used by Melling and Lewis (19S2) and observed values of currents 
adjacent to polynas reported by Schumacher et al. (19S3). The value of V 0 determines 
the areal flux used in calculating the salinity increase in the water column. The initial 
velocity is increased after the non-entraining, non-rotating rudiment plume travels 10 
km straight downslope. Effects of rotation are neglected at this early stage of the plume 
to minimize the exposure to brine rejection under the polyna. Increasing the velocity 
after 10 km ensures more realistic salinity values of the plume, and would, in any case, 
occur in the first steps of the streamtube model. The non-entraining, non-rotating 
rudiment plume continues straight downslope until a depth of 45 meters is reached. This 
depth was chosen due to its use by Melling and Lewis (19S2) as the initial depth of the 
streamtube model. 

At 45 meters, the plume descends the slope according to the dynamics set forth 
and established by Smith (1975), Killworth (1977). and Melling and Lewis (19S2). Three 
stratified layers are defined in this coupled model. Layer 1 lies between 45 and 55 meters 
depth. Layer 2 lies between 55 and 300 meters, and Layer 3 lies below 300 meters. The 
plume is assumed to interleave with ambient waters when the velocity is less than .01 
m's: however, all sensitivity studies will be carried out using 90 days as the maximum 
number of days the plume can descend while retaining its integrity. 

B. MODEL PARAMETERS 
1. Environmental variables 

In order to examine the environment's influence upon possible polyna- related 
production of halocline and deep waters, certain environmental variables were individ- 
ually varied during model simulations. The sensitivity of the model to a variable's al- 
teration is defined by the dependence of the plume's depth after 90 days upon the 
particular variable, relative to a standard case. The relationship between the depth of 
the plume and a single environmental variable will be reported in Chapter 3. The vari- 
ables to be investigated are listed in Table I on page 14. 
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Figure 2. Schematic representation of the coupled model. As the sheifwater is ex- 
posed to the polyna, a high density plume is formed, resulting in the 
plume's descent down the slope until it is neutrally buoyant. 
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Table 1 . LIST OF ENVIRONMENTAL INPUTS: The dependency of the plume 

depth after 90 days on each of the variables below will be examined in 
Section 3. 



Variable 


Definition. 


T.„ 


air temperature 


r. 


wind velocity 


s 


shelfwater salinity 


l\ 


initial shelfwater velocity 


Sw 


shelf width, from the 45 m isobath 




initial shelf slope 


0, 


secondary slope 


•V, 


Briint- Vaisala frequency, layer 1 


.Y, 


Briint-Vaisala frequency, layer 2 


AZ : 


thickness of layer 2 



2. Environmental conditions for the standard case 

To achieve consistency with pre-existing works, many of the environmetal con- 
ditions for the standard case were borrowed from Melling and Lewis (1982), Killworth 
and Smith (19S4). and Pease (19S7). The standard case initial values are listed in 
Table 2 on page 15. 

The most significant departure from an initial value given by previous work, is 
the H, value. In this paper, the frazil ice collection thickness is taken to be .2 meters. 
Pease (19S7) shows a higher collection thickness results in longer polyna opening time 
and larger polyna width for the same wind speed and air temperature. Weeks and 
Ackley (CRREL Monograph S2-1. 1 9 S 2 ) give higher values for frazil ice collection 
thicknesses in a wind driven scenario vice .01 to .1 meters in quiescent conditions. In 
addition to the frazil ice collection thickness, the salinity of the ice is set to 7 ppt (Cox 
and Weeks, 1974) instead of 5 ppt used in Killworth and Smith (1984). 
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Table 2. INITIAL VALUES USED FOR THE STANDARD CASE. 



Variable 


Value 


Definition 


T,„ "C 


-15. 


air temperature 


T,. "C 


-1.8 


water temperature 


I sec 1 


15. 


wind velocity 


f *. m sec -1 


3 ° o r,„ 


ice floe velocity 


II . m 


> 1 


frazil collection thickness 


o. 


5.67a* 10 _ - 


* Stephan- Boltzmann constant 


e. 


0.95 


emissivity of the air 


p c kgnr 2 


1.30 


air density 


p ? . kgnr 2 


1.026-vlO 3 


seawater density 


p.kgnr 3 


0.95xl0 3 


ice density 


c 


2.0-vlO- 3 


sensible heat coefficient 


C, Jdeg-'kg-' 


1004 


specific heat of air 


<X. iv»r 2 


301 


longwave radiation upward 


L. Jkg-' 


3.34x10 s 


latent heat of fusion 


b, . m 


10. 


flow thickness 


u/, ; . m 


1000. 


flow width 


1 m sec -1 


.04 


flow speed 


/? . degrees 


29 


flow direction 


/. sec- 1 


1.3Sxl0-* 


Coriolis parameter 


.V. sec -1 


0.0316 


Brunt- Vaisaia frequency, layer 1 


N 


0.0077 


Brunt- Vaisiila frequency, layer 2 


N 


0.001 


Brunt-Vaisala frequency, layer 3 


AZ, . meters 


195 


Layer 2 thickness 


0. 


0.5.rl O' 3 


initial bottom slope 


e 3 


5xl0- 3 


secondary bottom slope 


K 


0.01 


drag coefficient 


E 


0.072 sin 6 sin p 


entrainment coefficient 


S . ppt 


32.5 


initial shelfwater salinity 


s, ■ ppt 


7.0 


Salinity of frazil ice 
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III. RESULTS 



A. STANDARD CASE 

In order to examine the environmental influence on the possible polyna production 
of haloeline and deep waters, a standard case was defined with values which are typical 
of Arctic Ocean observations. The offshore wind velocity requires a considerable 
meteorological event providing winds higher than the monthly mean winds. However, 
the value of 15 m s for offshore wind velocity is within observed wind speeds for Arctic 
storms reported in previous papers (Schumacher, et al. (19S3) and Aagaard. et al. 
(1985)). The standard case plume obtains a salinity of 34.98 ppt after transversing the 
polyna. which corresponds to a density of 1032.93 gmjm 3 and a negative buoyancy of 
.0643 m/s 2 . It is therefore not suprising when the plume penetrates the same Arctic 
pycnocline used in Melling and Lewis (1982) to a depth of 436 meters. Figure 3 on page 
17 shows the velocity of the plume as it descends the slope. Due to the plume's high 
negative buoyancy, the plume accelerates slightly on the initial slope then rapidly in- 
creases in speed after the shelf break at 55 meters. The maximum velocity is 0.30 m s 
reached at a depth of 60 meters. The mass flux, AV, increased by over 1600% empha- 
sizing the dependence of the entrainment rate on the velocity of the plume. At the end 
of 90 days, the standard case plume has a velocity of 0.017 m s and a negative buoyancy 
of 0.000482 ml sec 2 . The plume has travelled 253 kilometers along the shoreline and 96 
kilometers down slope. The flow path of the plume can be seen in Figure 4 on page 
18. Overall, the plume in the standard case behaves similarly, although with different 
magnitudes, as the plumes in Melling and Lewis (1982). With this in mind, sensitivity 
studies of the impact of individual environmental parameters on the 90-day plume depth 
can be examined. Sensitivity experimental results typically include 10-20 data points with 
increased resolution in areas of abrubt change. 

B. SENSITIVITY EXPERIMENTS 

The polyna's maximum width is dependent on the growth rate of ice, the advection 
rate of ice, and the collection thickness of ice. Air temperature and the wind speed af- 
fects the growth rate of ice, while the advection rate is proportional to the wind speed. 
Lowering the air temperature results in a monotonic increase in the plume's depth at 90 
days as shown in Figure 5 on page 20. This is expected due to the increased production 
of frazil ice resulting in higher salinization of the water column. Although Pease (19S7) 



16 



LT) 




PLUME DEPTH (M) 



Figure 3. The standard case plume's velocity is shown with depth: Note the rapid 

increase in velocity as the plume crosses the shelibreak at the 55 meter 
isobath to a steeper slope. 

has shown the polyna size decreases with lowering air temperature at a given wind speed, 
the plume's salinity (hence, density) increases most rapidly in the first 10 kilometers of 
downslope travel. Here the plume will be subjected to significantly higher salinization 
due to higher ice growth rate and lower areal flux rate. Although the polyna size is 
smaller with lowering air temperatures, the majority of the salinization has already oc- 
curred in the first 10 kilometers. Indeed, the polyna shows a decreasing salt contribution 
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Figure 4. The standard case plume non path in the x-y coordinates: The effects 

of Coriolis can be seen as the plume's angle P decreases during the 
plume's descent. 



IS 



after the first 10 kilometers with decreasing air temperature. This can be seen in 
Table 3 on page 19. This decreasing salt contribution to the water column is due to the 
polyna's smaller width. 



Table 3. SALINITY CONTRIBUTION OF POLYNA WITH DECREASING AIR 
TEMPERATURE: The table shows the initial salinity of the plume. S,.„, 

after the first 10 kilometers of downslope travel, the final salinity of the 
plume. S„ after travelling the entire width of the plume. A'., and the de- 
creasing salt contribution of the polyna during the distance X p — \0km with 
decreasing air temperature. T .. r . - 



s„„ . ppt 


5, . ppt 


AS , ppt 


T ° C 

1 cir* 


33.36 


34.36 


1.0 


-5.0 


33.75 


34.67 


0.92 


-10. 


34.14 


34.9S 


0.S4 


-15. 


34.54 


35.30 


0.76 


-20 


34.93 


35.61 


0.6S 


-25 


35.32 


35.92 


0.60 


-30 


35.72 


36.23 


0.51 


-35 



On the other hand, if air temperature is held constant, and offshore wind speed is 
increased, the polyna grows in size due to increased ice advection rate. The initial con- 
tribution of salt during the first 10 kilometers of the polyna does not vary however, the 
increased polyna size exposes the water column to further salinization. This results in 
a denser, deeper penetrating plume with increasing wind speed. (See Figure 6 on page 
21). The dependence of plume depth upon the wind velocity and air temperature illus- 
trate the possible atmospheric influences upon the Arctic Ocean's salinity characteristics. 
The frequency and severity of Arctic storms which would be conducive to high density 
plume formation is a consideration in the calculation of production rates, but cannot 
be discussed quatitatively using this simplified model. 

The effect of varying the initial salinity of the shelfwater on plume depth is shown 
in Figure 7 on page 22. The results of this experiment show deeper penetration with 
higher initial salinity. This is not suprising due to the significant role of salinity upon 
water density near the freezing point. The plume depth after 90 days is greater than the 
depths reported by Melling and Lewis (1982). even with relatively low salinity values. 

Initial shelfwater velocity defines two parameters within the model, the areal flux 
used in computing the new salinity of the water column and the initial velocity of the 



19 



o 




Figure 5. Effect of decreasing air temperature on the plume's 90-day depth: With 

decreasing air temperature, the increased ice growth rate results in higher 
initial salinization of the water column. 

plume. The higher the initial velocity, the less saline the plume. However, when using 
the plume depth at 90 days as a standard criteria to examine the results, the higher ve- 
locity plumes reach deeper depths at 90 days titan the denser, but lower velocity plumes. 
If a speed criteria were used, then there would be a negative correlation between the 
plume depth and initial shelfwater velocity. Instead, Figure S on page 23 shows a pos- 
itive correlation between plume depth and initial shelfwater velocity. 



20 



o 




Figure 6. Effect of varying offshore Mind speed on the 90-day plume depth: The 

polyna's increased width with increasing wind speed results in a higher 
density plume due to longer exposure of the water column to the 
polyna's brine rejection. This results in a deeper penetrating plume. 

Similarly, the parameters of shelf slope and secondary slope appear to have contra- 
dictory results. Figure 9 on page 24 shows as the initial shelf slope is lessened, the 
plume penetrates to deeper depths after 90 days of travel. In this experiment, the shelf 
width was fixed at 20 kilometers after the 45 meter isobath. This deeper penetration 
appears to be due to a lower entrainment parameter, which is a function of the shelf 
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Figure 7. The effect of initial shelfwater salinity on the 90-day plume depth: With 

increasing shelfwater salinity the plume attains a higher density due to 
increased salt rejected by ice growth added to the shelfwater. 

slope. In contrast. Figure 10 on page 25 shows deeper penetration with an increased 
secondary slope. Although the entrainment parameter still increases with slope, the 
90-day plume depth is limited by the slope itself Since 90 days was considered a rea- 
sonable time scale for the plume to maintain its integrity, deeper plume penetrations 
would most likely occur with steeper slopes after the shelf break, i.e. in submarine 
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Figure 8. Effect of varying initial shelfuater velocity on 90-day plume depth: An 

increased initial velocity results in a deeper penetrating plume at 90 days. 

canyons. Varying the shelf width after the 45 meter isobath has little effect on the 
90-day plume depth. This is shown in Figure 1 1 on page 26. 

Examining the stratification of the Arctic Ocean is accomplished by varying the 
Brunt- Yiiisala frequency of Layers 1 and 2. The deep Layer 3 remains constant at a 
value used in Killworth (1977) for deep waters. Layer 1 is the initial pycnoc.line and re- 
presents the transition from the fresh surface layer to the saltier halocline. Figure 12 
on page 27 shows the 90 day plume depth is little changed even with a substantially 
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Figure 9. This figure shows the effect of the initial shelf slope on the 90-day plume 
depth: Due to the entrainment parameter being a function of slope, a 

steeper slope results in higher entrainment of ambient water, thus more 
rapid loss in negative buoyancy. 

weakened initial stratification with A" 2 = 2.0.tl0~ 4 sec 2 . The short length of time and 
space which the plume is exposed to this layer is probably the cause of its diminished 
impact on the 90-day plume depth. 
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Figure 10. The effect of varying the secondary slope on the 90-day plume 
depth: Even with a higher entrainment rate, the plume reaches a 

deeper depth with a steeper slope, due to the physical boundary of the 
ocean bottom. 

On the contrary. Figure 13 on page 28 shows the dramatic difference in the 90-day 
plume depth with small changes in the Brunt- Vaisala frequency in this 2-15 meter thick 
layer. Additionally, if the second layer's thickness is decreased, the plume will penetrate 
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Figure 11. The effect of varying the shelf width (after the 45 m isobath) on the 
90-day plume depth: Varying the shelf width has little effect on the 

90-day plume depth. 

to deeper depths. This is illustrated in Figure 14 on page 29. Note the discontinuity in 
the slope of the curve is caused when the plume no longer pentrates the second layer. 
Next, the standard case is compared in detail with a "submarine canyon" simulation. 
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Figure 12. The effect of varying the initial pycnociine strength on the 90-day plume 
depth: The initial pycnociine. between 45 and 55 meter depths, has a 

small elTect on the 90-day plume depth due to the small thickness of 
this layer. 

C. THE "SUBMARINE CANYON" CASE 

The "submarine canyon" case is a modification of the standard case using a steeper 
secondary slope, lower Layer 2 Briint-Vaisalii frequency value, and a decreased Layer 2 
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Figure 13. The effect of varying Layer 2 Brunt- Yiiisala frequency on the 90-day 
plume depth: Due to the plume's longer exposure time in this layer, 

changes in Briint-Vaisala frequency tend to have a more dramatic effect 
on the vertical motion of the plume. 

thickness. In addition, the angle /?. is forced to the initial value of 29 degrees to simulate 
the steering of the plume by topography. The initial values used in the "submarine 
canyon" case are shown in Table 5 on page 31. Forcing the angle /? towards its initial 
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thickness, the 90-day plume depth becomes constant. 



value results in the plume travelling down the V axis further and less down the X axis. 
Of course, the value chosen for forcing /? will alter the X-Y flow path of the plume. The 
How path of the "submarine canyon" case is seen in Figure 15 on page 32. The plume, 
descends to a depth of over 1300 meters in 90 days, has moved 207 kilometers along- 
shore and 160 kilometers off-shore. In Figure 16 on page 33. the "submarine canyon" 
plume is seen to reach a maximum velocity of 0.36 m s whereas the standard case plume 
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has a maximum velocity of 0.30 m s (Figure 3 on page 17). The effect of increased 
slope, the angle /?, and velocity on entrainment is seen when comparing the volume flux 
of the "submarine canyon" plume (Figure 17 on page 34) and the standard case plume 
(Figure 18 on page 35). 

D. RESULTS USING OBSERVED TEMPERATURE AND SALINITY FIELDS 

The values for this experiment are the same as the standard case except the Briint- 
Vaisalii frequency is computed using observed values of temperature and salinity. The 
temperature and salinity records are from Ostlund (19S7) and Aagaard et al. (1981) 
which were taken in the Canadian Basin. The results of the experiment are shown in 
Table 4. The 90-day plume depth and maximum velocities are very similar to the results 
from the Layer 2 Briint-Vaisala frequency sensitivity experiment. 

In an attempt to model a current as observed in Aagaard et al. (1985). a plume 
width of 20 kilometers was used together with a steep secondary 5 slope and the 
temperature salinity fields reported in Aagaard et al. (1981). The plume's velocity with 
increasing depth is shown in Figure 19 on page 36. Aagaard et al. (1985) reported a 25 
kilometer wide. 15 meter deep plume with a mean velocity of 0.45 m s. The simulation 
yields a velocity of 0.30 m s before the shelf break and 0.92 m, s while descending the 
steep secondary slope. 



Table 4. THIS TABLE COMPARES SOME RESULTS FROM USING ACTUAL 
T-S RECORDS: The AIWEX-2 and AIWEX-8 stations are from 

Ostlund( 19S7) and the third record is taken from Aagaard et al. (19S1). 



Variable 


AIWEX-2 


AIWEX-8 


Aagaard 


^ max m S 


.397 


.3S9 


.395 


90-day depth, m 


217 


213 


305 


X distance, km 


311 


312 


363 


'll' distance, km 


380 


375 


478 
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Table 5. INITIAL VALUES USED FOR THE "SUBMARINE CANYON" 
CASE: The highlighted values are those parameters which differ from the 

standard case. 



Variable 


Value 


Definition 


T, ... C 


-15. 


air temperature 


T„ . "C 


-l.S 


water temperature 


tn sec* 1 


15. 


wind velocity 


V. m sec* 1 


3 % i;„ 


, ice tloe velocity 


a . m 


*> 


frazil collection thickness 


a, lVnr : deg~ 4 


5.67x10-* 


Stephan- Boltzmann constant 




0.95 


emissivity of the air 


p 7 kgm~ 2 


1.30 


air density 


kg nr* 


1.026.vl0 3 


seawater density 


p.kgnr* 


0.95.V 10 3 


ice density 


C, 


2.0-rlO- 3 


sensible heat coefficient 


C. Jdcg-'kg- 1 


1004 


specific heat of air 


Q: u . IVnr 2 


301 


longwave radiation upward 


L. Jkg 


3.34-vlO 5 


latent heat of fusion 


/?., . m 


10. 


flow thickness 


;r, , m 


1000. 


How width 


( m sec* 1 


.04 


How speed 


ft . degrees 


29 


flow direction 


/. sec -1 


1.33.rlO- J 


Coriolis parameter 


.V. sec -1 


0.0316 


Brunt-Vaisala frequency, layer 1 


N 


0.006236 


Brunt-Vaisala frequency, layer 2 


A Z, . meters 


145 


Layer 2 thickness 


X 


0.001 


Brunt-Vaisala frequency, layer 3 


0, 


0.5.Y10- 3 


initial bottom slope 


0, 


9y10- 3 


secondary bottom slope 


K 


0.01 


drag coefficient 


E 


0.072 sin 6 sin (I 


entrainment coefficient 


5 . ppt 


32.5 


initial shelfwater salinity 


5, . ppt 


7.0 


Salinity of frazil ice 
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forced to the initial /> value to simulate the topographic steering within 



a submarine canyon. 
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Figure 16 . The velocity of a "submarine canyon" plume with increasing depth: The 

plume has the characteristic "jump" in velocity after the shelf break 
where it reaches a maximum velocity of 0.36 m. s. 



33 



o 




Figure 17. The normalized mass flux of the "submarine canyon" plume: The 

plume increases in mass by over 6S00% indicating only 1 5% of the 
water reaching the deeper depths is shelf water. 



34 




35 



r\j 




PLUME DEPTH ( M ) 



Figure 19. Velocity shown with depth of a 20 kilometer wide plume: This figure 

shows a wider plume can penetrate deeper than the narrower 1 kilo- 
meter plumes used in the previous experiments. 
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IV. DISCUSSION 



This paper has attempted to provide a mechanism in which high salinity, low tem- 
perature water is produced and transported otT the continental shelves into the Arctic 
halocline or deeper depths. The environmental parameters which affect the production 
and transport of this cold, salty water were examined, resulting in the conclusion that a 
polyna can produce water with a salinity greater than 34.5 ppt near the freezing point 
which can be transported via gravity flow plumes. The results show deeper penetrations 
to depths greater than 700 meters is possible depending on the set of environmental pa- 
rameters. This indicates polyna-produced high density plumes may also be responsible 
for the ventilation of the deep Arctic Basins, specifically the Canadian Basin. 

Ostlund et al. { 1 9S7) placed the residence time of Canadian Basin Deep Water be- 
tween 500 and 800 years. Using the mixing ratio of 2:1 (shelfwater to intermediate wa- 
ter) from Aagaard et al. (19S5), shelfwater with a salinity of 35.1 ppt requires an 
off-shelf flux of approximately .0063 Sv to ventilate the Canadian Basin in 800 years. 
The results from the 20 kilometer wide plume experiment give an annualized flux of 
0.0035 Sv. In other words, an average of about two such events per year where 
polyna-produced plumes penetrate below the temperature maximum over the past 800 
years could account for the ventilation of the Canadian Basin. Furthermore, Ostlund 
et al. (1987) states that the deep water probably consists of much less than 10% 
shelfwater and would have to be a salinity of 36.15 ppt. The "submarine canyon" case 
shows the plume reaching a depth of over 1300 meters consisting of only 1.3% 
shelfwater at this depth. Table 3 on page 19 shows under extremely cold conditions, a 
salinity value of 36.23 ppt can be achieved with this model. Therefore, it is reasonable 
to expect that polyna-produced saline water may significantly contribute to the venti- 
lation of the Canadian Basin. 

Climatological variations which would affect the number of deep penetrating plumes 
are worthy of discussion. Besides the major periods of global ice ages and their 
counterpart warm eras, there have been recorded several shorter cold and warm epochs. 
These shorter periods coincide with weak and strong wind circulation, respectively. In 
the 20th century, both a warm period lasting about 20 years and a colder period begin- 
ning around 1950 have been recorded. In each of these two eras, there has been a sud- 
den outflow of surface layers from the Arctic Ocean. Increased ice flow was recorded 
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in 1938 around Iceland during the warm period, while an increased flow of cold, low 
salinity water reached Iceland from the Arctic in the 1%0's. Removal of the colder, 
fresher surface layer from the Arctic Basin would reduce the vertical temperature and 
salinity gradients resulting in a less stable vertical structure for the ocean. In effect, an 
environment similar to reducing the model's Layer 2 Brunt-Viiisala frequency and re- 
ducing the Layer 2 thickness may occur, enhancing deep plume penetrations. In addi- 
tion to possible weakening of the vertical structure, removal of the ice cover leaves large 
areas of open water where heat is more readily removed from the ocean. In these large 
areas of open water, frazil ice grows rapidly and abundantly to increase the brine rejected 
into the water column. In the warm period of the early 20th century, the extent of the 
Arctic's semi-permanent ice pack was greatly reduced (up to 10%) between 1920 and 
1938. The reduction in the semi-permanent ice pack is attributed to higher ice removal 
rates (increased advection) due to the strong wind circulation of the period. In Ostlund 
et al. (1987). it is noted that the deepest waters in the Canadian Basin may be the re- 
mains of the warm climate which occurred around 1000 A.D. vice any cold epochs such 
as the "Little Ice Age” in the late 17th and early 18th centuries. 

In summary, periods when Arctic wind circulation is strengthened resulting in open 
water due to rapid ice divergence coupled with a weaker vertical Arctic Ocean stability 
seems to be the most conducive environmental scenario for deep plume penetrations to 
provide the additional salt source in the Canadian Basin. 

In Killworth and Smith (19S4), the heat diffusivity constant, k , required to achieve 
steady state balances requires some peculiar changes in magnitude with increasing depth. 
A value of k^3 x 10 -4 m\s _1 is established to achieve a steady balance between the non- 
turbulent shelfwater plume used in Killworth and Smith (1984) and the region where the 
plume interleaves (note: the density of this plume is based on the brine rejection due to 
a single year's ice growth over a continental shelf of 50 meters average depth). To 
properly supply the East Greenland current outflow, k was required to be 
so x 10 _ -‘m 2 sec -1 . And finally, to balance the inflow of warm Atlantic water and upward 
advection of heat, k at this level was required to be ^1.2 x 10' 4 m 2 sec* 1 If the normalized 
mass flux of the turbulent plume model was to be used as a gauge to understand the 
turbulent diffusion, this thesis would present some interesting results. In Melling and 
Lewis (1982), the normalized mass flux was of order ^10°. In the standard case, the 
normalized mass flux was of order ^10 l . The "submarine canyon" case had a mass flux 
of order ^10 2 . Could it be that the three above cases provide a insight into the missing 
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physics so stated in Kilhvorth and Smith ( 1 9S4 >? Some other missing physics may be the 
double diffusion process as outlined in Carmack and Aagaard (1973) and the breaking 
of internal waves (Perkin and Lewis. 1978). Aagaard et al. (1985) suggests the incorpo- 
ration of turbulent plume dynamics into the "filling box" model of Kilhvorth and Smith 
(19S4). A mechanism providing a specific low temperature, high salinity water source 
has been developed in this paper. Although the mechanism is episodic, as suggested 
Aagaard et al. (19S5), (Kilhvorth and Smith (1984) had earlier used the term 
"spasmodic"), due to the dependence on environmental variables, flux rates of this highly 
saline water can be estimated from renewal times of the Canadian Basin. The problem 
of incorporating the turbulent plume in the "filling box" model may be tackled, but is 
beyond the scope of this paper. 

The maintenance of the Arctic halocline by the mechanism presented in this study 
alone would require many times the number of active polynas producing saline 
shelfwater. Using the 3:2 mixing ratio proposed by Aagaard et al. (1981). a minimum 
of 2.5 Sv production rate of cold, saline shelfwater was determined necessary to replenish 
the halocline over a period of 10 years. With the flux rate of the 20 km wide plume (.0035 
Sv), the individual mechanism modelled in this paper would require 200-300 coastal 
polynas producing 2-3 high salinity, low temperature plumes per year to maintain the 
Arctic halocline. Although 200-300 such polynas may be thought prohibitive, areas of 
open water which exist during the onset of freezing temperatures in the late fall, as well 
as leads which continue through winter over as much as 10% of the ice pack, could 
provide a rapid source of salt to help overcome the summertime salinity gradient. In 
times of strong vertical structure, these polyna-produced plumes can transport low- 
temperature. high salinity water into the Arctic halocline. Overall, the polyna seems to 
be an important source of brine which may help maintain the Arctic halocline and ven- 
tilate the deep basins. 
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